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Abstract. We study the harmonic analysis of the quadrature mirror filters 
coming from multiresolution wavelet analysis of compactly supported wavelets. 
It is known that those of these wavelets that come from third order polyno- 
mials are parametrized by the circle, and we compute that the corresponding 
filters generate irreducible mutually disjoint representations of of the Cuntz 
algebra O2 except at two points on the circle. One of the two exceptional 
points corresponds to the Haar wavelet and the other is the unique point on 
the circle where the father function defines a tight frame which is not an or- 
thonormal basis. At these two points the representation decomposes into two 
and three mutually disjoint irreducible representations, respectively, and the 
two representations at the Haar point are each unitarily equivalent to one of 
the three representations at the other singular point. 



1. Introduction 

In this paper we show that wavelets may be constructed from representations of 
two systems of operator relations, one on L 2 (K.) and one on L 2 (T), for the case of 
one real dimension. Focusing on the case of compact support, the analysis reduces 
to a certain finite-dimensional matrix problem which is especially amenable to an 
algorithmic and computational approach. The associated algorithms are worked 
out in detail for a variety of examples which includes the Daubechies wavelet, and 
which also reveals some perhaps unexpected symmetries. 

One benefit from the representation theoretic approach to wavelets is that it 
provides a coordinate-free way of making precise notions of irreducibility which oc- 
cur in the wavelet literature without always having precise definitions. Specifically, 
examples in L 2 (R rf ), for d > 1, may occasionally be reduced to simpler examples 
in one dimension, i.e., in L 2 (R), by a tensor product construction, but this analysis 
depends on the chosen spatial coordinates in M d , while the representation-theoretic 
approach in the present paper does not. 



One of our results, Corollary |3^, specifies in a general context (for compactly 
supported wavelets in M.") a decomposition formula (finite orthogonal sums of ir- 
reducible representations) for the representation associated with a system of high- 
pass/low-pass filters which generate the wavelets in question. 

It has been known for some time that a class of convolution operators from sig- 
nal analysis, called subband filters, satisfy certain operator relations |3l], Lemma 
2.1]. Perhaps it is less well known among experts in multiresolution wavelet theory 
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that these operator relations were introduced in C* -algebra theory by J. Dixmier 
|l4| Exemple 2.1] and J. Cuntz several decades ago, and the C*-algebra they 
generate is now called the Cuntz algebra of order N and is denoted by Oat, where 
N is the scale of the resolution. This algebra is independent of the particular scale- 
N multiresolution wavelet, but the unitary equivalence class of the corresponding 
representation may depend on the wavelet. The detailed structure of these rep- 
resentations has, however, so far only been worked out in the single case of the 
Haar wavelet (see below). The purpose of the present paper is to work out the 
structure of these representations for all compactly supported wavelets, using a 
method tailor-made for the purpose in M. We will show that all representations 
obtained from compactly supported wavelets have a finite-dimensional commutant, 
and as a consequence they decompose into a finite direct sum of irreducible repre- 
sentations. We also display a one-parameter family (with two singular points) of 
mutually inequivalent representations of O2 on L 2 (T) for which the correspond- 
ing family of wavelets contains Daubechies's continuous, one-sided differentiable 
mother function, ip g L 2 (R), supported on [0,3] C R. In our one-parameter fam- 
ily of wavelets supported on [0,3], there is actually a left-handed and a paired 
right-handed Daubechics wavelet, resulting from a natural symmetry in the fam- 
ily. In going from one to the other, the one-sided differentiability property reverses 
direction. 

Let us briefly review how one constructs representations from a multiresolution 
wavelet of scale N. Many more details may be found in [El. Excellent accounts of 
multiresolution wavelet analysis in general may be found in pl| and || . 

Define scaling by N on L 2 (R) as the unitary operator U given by (U£) (x) = 
N~i£ (7V _1 x) for ( e L 2 (R) , x G R, and translation as the unitary operator T 
given by (T£) (x) = £ (x — 1). There is a father function or scaling function ip which 
is a vector in L 2 (R) such that 

(1.1) {T k (p) k& is an orthonormal set in L 2 (R) . 
Furthermore, one assumes that there is a sequence (b n ) € i 2 such that 

(1.2) U<p = ^2b n T n (p, 

n 

and then necessarily \b n \ 2 — 1. (It seems to be fairly conventional in wavelet 
theory to only consider real &, but this is not too important for what follows.) 



A weaker, so-called "tight frame", propert y for the vectors in (1.1) will also be 



L 2 (R) spanned by {T k ip} k one also assumes 



considered as a degenerate case in Section [LL2| . If V is the closed subspace of 



(1.3) A unv « = w ' V unv ° = l2 ( r ) • 

These are all the properties of the father function cp that are needed. One 
example is the Haar father function ip (x) = X[o,i] ( x )- 
Define a function mo in L 2 (T) by 



(1.4) m (t) = m Q (e- u ) = ]T fo r 



e -int 
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Choose functions mi, . . . , thn-x in L 2 (T) such that 
N-l — 



(1.5) 



^ rrij I £ + 



fc=0 



2?rfc 



m, £ + 



27rfc 



for almost all £ G K, i, j = 0, 1, . . . , N — 1, or, equivalcntly, such that the N x N 
matrix 

/ m (z) mo(pz) ... m (p N ~ 1 z) \ 



(1.6) 



mi(z) mi(pz) 



\m N -i(z) rriN-i(pz) 



N-l, 



rriN-iifi Z )J 



where p = e^' , is unitary for almost all z G T. (With mo given as above, 
mi, . . . , iriN-i may always be so chosen; see, e.g., ||.) If we define . . . , ipN-i € 
i 2 (R) by 

(1.7) VN$ i (Nt) = m i (t)$(t) 

for £ G R, i = 1, ... ,7V — 1, where * denotes Fourier transform, unitarity of the 
above matrix is equivalent to orthonormality in L 2 (M) of the set 

{ C/n ^}„, fceZ;i = 1 ,...,iv-i • 
The V'j's are called the mother functions. If TV = 2, there is only one, of course. 

Unitarity of (1.6) is also equivalent to saying that the operators Si, defined on 
L 2 (T) by 

(1.9) (S i 0(z) = m i (z)t{z N ) 

for £ G L 2 (T), z G T, i = 0, 1, . . . , N - 1, satisfy the relations 



N-l 



(1.10) 



which are exactly the Cuntz relations. There is a one-to-one correspondence be- 
tween operator solutions to (1.1C) and representations of Oat, and since On is 
simple, these representations are always faithful. The Fourier transform of S* (the 
adjoint of ([1.9])), acting on i 2 (Z), is the quadrature mirror filter F{ in |3^]: Fq is 
low-pass, and Fx,... ,Fn-i are the corresponding high-pass filters for the signal 
reconstitution process. Let 



(1.11) 



(z)=5>«z" 



be the Fourier decomposition. It follows from (1.9) that for x = (xk) keZ G £ 2 , we 
have 



(1.12) 



U) 



n-Nk X ki 



( F 3 X )n =Ys a k-Nn X k, 

as operators t 2 — > t 2 . The Cuntz relations in £ 2 -operator form, 

N-l 

(1.13) ^=^1, ^FJF^l, 

3=0 

then summarize subband filtering, which can be written in diagram form as in 
Figure |l|. Here "analysis" is splitting into subbands and the application of Ft , and 
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Figure 1 . Signal subband filtering 



"synthesis" is the application of F* followed by summing over the subbands again. 
The low-pass subband corresponds to i = 0, and the high-pass subbands correspond 
to i = 1, . . . , N - 1. See |l) and § for details. 

The Oj\r-representations given in (l.S) play a crucial role in the wavelet analysis 



in a second related way. A scale- TV wavelet in L 2 (K) is an orthonormal basis (or a 
tight frame) of the form (l.S) as d escribed above. An important point is that the 
corresponding Si-operators of (|l.9| ), which constitute the ©Ar-representation, enter 
directly and explicitly into a formula for the L 2 (R)-cxpansion coefficients c n ki of 
£ = En,M °nki (0 U n T k ip h £ e i 2 (R), and we refer to |, eq. (1.35)] for details on 
that. 

that the scaled vectors [/-0i and C/y are both finite 
if and only if the functions rrn are 



We see from (1.7) and (fl 

linear combinations of translates {T fe <p} fcgi 
polynomials, and this is reflected in the fact that the wavelets </?, ipi have compact 
support if and only if all the functions rrii (z) are polynomials in z. (See p2[ Chapter 
5], |2l], Section 3.3].) In Q, a detailed study was made of the representations of 
On defined by fll.SP in the case where m, (z) are monomials (or more precisely, 
monomials of the form m, (z) = z ni \ the more general case where m; (z) — XiZ Ui 
with Xi e T C C was considered in |L3)). It is clear from (1.2) and (1.7) that 



the other mi-functions coming from wavelets are never monomials, bu t the Haar 
wavelet (for N = 2), <p(x) = X[o,i] ( x )> ^ s close: one checks from (1.2) and (1.4) 
that mo (z) = (1 + z) fy/2. The most general choice of mi is then 



(1.14) 



mi (z) = zf (z 2 ) m (-z), 



where / maps T into T, and one conventional choice is / = — 1, i.e., 
(1.15) mi (z) = (1 - z) /V2. 

Thus the Haar mother function is given by -75^ (f ) = (</? (%) — ¥ ( x — 1)), i-e. , 
the graph of tp is that represented in Figure 
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Figure 2. Father and mother functions for the Haar wavelet 



If Si is defined by fll.9|), and one transforms the representation by A- ( \ ^ ) S 
U(2), i.e., 

(1.16) T = (5o + S!)/A T^tSo-SjJ/VS, 

one verifies that the pair To, T\ still satisfies the Cuntz relations, and 

(1.17) T Q Z(z)=t(z 2 ), T^(z) = z^(z 2 ). 

This is one of the monomial representations studied in H, and by |^|, Proposi- 
tion 8.1], this representation of O2 decomposes into two inequivalent irreducible 
subrepresentations on the subspaces 

(1.18) H 2 (T) = span 11 ' ^ {z 11 \ n € N U {0}} , 

(1.19) H 2 (T) 1 - = zH 2 (T) = spall 11 ' "2 [ z ~ n \ n e N} , 

where zi/ 2 (T) refers to complex conjugation, and " ' " 2 is L 2 (T) closure. Thus 
the original Haar wavelet representation is a direct sum of two inequivalent irre- 
ducible subrepresentations. In general, when the functions rrii are polynomials, this 
simple trick of reducing to monomials is not going to work, but we will see that 
it is possible to develop a theory for polynomial representations which nonetheless 
has many general analogues with the monomial theory. 

2. Finitely correlated states on the Cuntz algebra On 

Let us recall a few facts about the Cuntz algebra O n from (To) , and the part of 
the results from || that will be needed in the sequel. 

If N G {2,3, .. . }, the Cuntz algebra On is the universal C*-algebra generated 
by elements sq, . . . , sn-i subject to the relations 

(2.1) s*Sj = 6ijl, ^2 s j s *j = 1 - 

The Cuntz algebra may be viewed as an interpolation between the algebra of 
the canonical anti-commutation relations (CAR) and the algebra of the canonical 
commutation relations (CCR): The g-canonical commutation relations, 

fljfl* - qa*a l = 5^1, 

i,j = 1, . . .,d, reduce to the CCR relations if q = 1, the CAR relations if q = — 1, 
and the Cuntz relations ( fO| ) if q = 0. See 0, |, [l7|] for details on this. 

The Cuntz algebra is a simple separable C* -algebra not isomorphic to the algebra 
of compact operators on a Hilbert space. Therefore the space of unitary equivalence 
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classes of irreducible representations of On cannot be parametrized in a measurable 
way |l5| . In this paper we will show that the representations coming from low-pass 
filters of genus 2 form a (necessarily tiny) one-dimensional variety in this enormous 
space. 

There is a canonical action of the group U (N) of unitary N x N matrices on 
On given by 

(2.2) T g (si) = W* s j 

for g = [gij] G XJ(N). In particular the gauge action is defined by t z (si) = zsi, 
z G T C C . If UHFjv is the fixed point subalgebra under the gauge action, then 
UHFat is the closure of the linear span of all Wick ordered monomials of the form 
s ii ' " S ik S jk" ' s ji- UHFjy is isomorphic to the UHF-algebra of Glimm type N°°, 

oo 

(2.3) UHFjv = Mn°° = <S> M n , 

l 

in such a way that the isomorphism carries the aforementioned Wick ordered mono- 
mial, Si ■ ■ ■ Si s* ■ ■ ■ s* , into the matrix element 

(2.4) e iin (g) e i3 j 2 (g) • • • <g) e ik3k (g) 1 <g> 1 (g) • • • . 
The restriction of r g to UHFjy is then carried into the action 

(2.5) Ad (g) ® Ad (g) ® • • • 

on <^)^° Mn- We define the canonical endomorphism A on UHFat (or on On) by 

(2.6) \(x)=J2 S 3 XS *3 

jez N 

and the isomorphism carries A over into the one-sided shift 

(2.7) Xi ® X2 ® X3 (g • • • > 1 (g Xi (g X2 Ig • • • 

on 0fMjv. 

If Si i— * Si € B (H) is a representation of the Cuntz relations on a Hilbert space Tt, 
we will say (by abuse of terminology) that the representation is finitely correlated 
if there exists a finite-dimensional subspace /C C H. with the two properties 

(2.8) S*)C C JC for i £ Zjv, 

(2.9) K, is cyclic for the representation Si i — > Si. 

The presence of such a finite-dimensional subspace /C is a special property of each 
of the representations under discussion, and therefore of the states of Oat which 
correspond to the representations. These states were studied in || with a view to 
the present applications. 

If P: Ti — * K. is the orthogonal projection onto /C, then ( |2.8| ) can be formulated 

as 

(2.10) Vi = PSi = PSiP. 
If we view Vi as operators in B (JC), we have 

(2.11) Y, V i V * = 1 - 
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and conversely, if Vi are operators in B (K) satisfying ( 2.1l[), they determine a 
representation Sj i— > Si of the Cuntz relations such that ( 2.1 0| ) is valid, and this 
representation is unique up to unitary equivalence if we require K to be cyclic |j| 
Theorem 5.1]. 

If JC\ is another Hilbert space and Wo, . . . , Wn-i are operators on Ki satisfying 

(2.12) W i W *i = a > 

and Si i— > Ti is the associated representation of On, then there is an isometric linear 
isomorphism between intertwiners U : Tiv — * Ti-w, i.e., operators satisfying 

(2.13) USi = TiU, 
and operators V G B(JC,JCi) such that 

(2.14) p(V) = WiVV* = V. 

This linear isomorphism is given by 

(2.15) U i — ► V = PiVP, 

where Pi : Hw ~ ¥ £i is the orthogonal projection onto K\. All these results do not 
depend on K and K\ being finite-dimensional, and they are given in || Theorem 
5.1]. 

An important special case is K\ = K and W% = Vj . Then p is a completely 
positive unital map, and the linear isomorphism ( 2.15| ) is an order isomorphism 
between the fixed point set of p (which is not necessarily an algebra) and the 
commutant {5^, S* \ i €E T^n}' ■ In particular, we have the following principle. 

(2.16) The representation s$ i — ► Si is irreducible if and only if p is ergodic: 
{AeB(fC)\p(A) = A} = CI. 

The rest of the discussion in this section can only be partially extended to the case 
when fC is infinite-dimensional (see f| S ection 6] for details). Define cr = p in the 
case when K,\ = K, and Wi — Vi in ( 2.14 ). If a is ergodic, then B (fC) has a unique 
cr-invariant state w. This state need not be faithful (see the example after the proof 
of Lemma 3.4 in ||). If E is the support projection of <p, then S*EK. C EK, for all 
i G 1>n (see ^ Lemma 6.1]). In that case, replace P by E, Vi by EVi, <x by the <r 
defined by the new Vi's on EK,, and then define a state ip on On by 

(2.17) <l>(S I S}) = <p(ES I S*jE). 

It was proved in Theorem 6.3] that the following three subsets of the circle group 
T are equal: 

(2.18) {t G T | ip o n = ip}, where r is the gauge action; 

(2.19) {t G T | t/j o Tf is quasi-equivalent to ip}; 

(2.20) PSp (cr) n T, where PSp (cr) is the set of eigenvalues of cr. 

(Of course, in the present setting, where EK. is finite-dimensional, PSp (cr) = 
Sp (cr).) Furthermore, this subset is a finite subgroup of T. If k is the order of 
this subgroup, the restriction of the representation to UHFjv decomposes into k 
mutually disjoint irreducible representations, and these are mapped cyclically one 
into another by the one-sided shift A. More specifically, one has PSp (cr) n T = 
PSp (A) n T, and, if tk = e~ , there exists a unitary U on TL, unique up to a scalar, 
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implementing Tt k , and such that U = 1. The operator U is the unique (up to a 
scalar) eigen-element such that A (U) = ikU. If 

(2-21) U =J2 f k E n 

is the spectral decomposition of U, then the spectral projections E n project into 
mutually disjoint irreducible subspaces invariant for the representation restricted 
to UHFjv, and A (E n ) = E n+ i, with A extended to B (Ti) by the formula A ( • ) = 



3. Polynomial representations 



From the relation (l.S) it follows that 

(3.i) (st$(z) = ± J2 ^M£M, 

W N =Z 

where the sum ranges over all JV'th roots w of z || eq. (1.17)]. Recall that the 
Fourier series version of (3.1) on £ 2 (Z) is the filter operator Fi of (1.12). In order 
to incorporate the monomial results obtained in 0, and also to make the present 
results applicable to wavelets in dimension v > 1, let us extend the definitions of the 
representations somewhat. We replace L 2 (T) with L 2 (T v ) and fix a matrix N with 
integer coefficients such that |det (N)| = N € {2, 3, . . . }. If z = {z x , . . . , z v ) € TT 
define 



(3.2) 



,N 



> *1 



if N = \nij\"j_ v (Note that this definition of z N is different from the one after 



(1.8) in B. The present convention implies that relations like [z 



N 



M 



= z NM and 
below. The 

is the transpose of the map x i— ► ~Nx on W passed to the 



are valid, where z n is defined as in connection with rt3 



(z N )" 

present map z 

quotient TT = W /2irW .) The map z i-> z N is AT-to-1. Let <t , . . .,ojv-i denote 
sections of this map, i.e., each m : T v — ► T" is injective, /i (ct, (T") n <7j (T w )) = 
if i 7^ j, where fi is normalized Haar measure on T v , and /i (<7j (Y)) = j^fi (Y) for 
all Borel sets FcT". Thus Ui<?z Gi C^") = up to sets of measure zero. The 



unitarity condition (1.6) then says that the N x N matrix 



(3.3) -= 



( m Q (<To(z)) 
m 1 (a (z)) 



m ((Ti(z)) 
mi(cri(z)) 



\m N -i(ao(z)) m N -i(ai(z)) 



m (a N -i{z)) ^ 
mi{a N -\{z)) 



rriN-i(crN-i(z)) J 



is unitary for almost all z £ T u . The representation (1.9), (3.1) of On now takes 
the form 



(3.4) 

and then 



{SiO (z) = mi (z)ti(z N ), 



(3.5) 



(S*0(z) = - ^ rrn (w)£ (w) 



COMPACTLY SUPPORTED WAVELETS AND CUNTZ RELATIONS 



9 



Now, assume in addition to unitarity of (3.3) that mo, . . . , m^-i all are polyno- 
mials, so that there exists a fixed finite subset DcZ" such that 

(3.6) m 3 (z) = ^z n . 

neD 



Here we have used the notation z n — . . . , z v ) 



(n,i,...,n„) _ rii n 2 
— z ± z 2 



2™" , and 



a„ e C. Let e„, n G Z u , denote the usual Fourier basis for L 2 (Z"), i.e., e„ (z) 
z™. It follows from (3^) that 



(3.7) 



Si). 



keD 



If in general we define ai' — when k ^ D, it follows from (3.7) or (3.5) that 



(3.8) 



*K = E < 



0) 



in 

Thus both Sj and 5* map trigonometric polynomials into trigonometric polyno- 
mials in this case. If the matrix N _1 defines a contractive map — >■ M" in some 
norm, one can say more. The following proposition is an analogue of Lemma 3.8 
in Q in the present setting. 

Proposition 3.1. Assume that all the (complex) eigenvalues of N have modulus 
greater than 1. It follows that there is a finite subset H C V with the property that 
for any n € 1 V there exists onMeN such that 

(3.9) S}e n £ f[H) = span {e m \ m e H} 
for all multi-indices I with \I\ > M. 

Proof. Let us give two proofs of this statement, both based on a study of the maps 
o~p : M. u — > W defined for p e D by 

(3.10) a p (x) =N- 1 (x-p) 

for By considering a Jordan form of N, as in the proof of Lemma 3.8 in Q, 

the condition | A,; | > 1 on the eigenvalues of N means that there exists a norm on C" 
such that ||N _1 |< 1 in the associated norm on B (C v ). If d = max{||p|| | p 6 D} 7 
it follows from ( 3.10 ) that \\a p (x)\\ < ||N _1 || (||x|| + d) for p € D, and by iteration, 

n 

(3.11) \\o- pi a P2 ---a Pn (x)\\ < \\x\\ + E || N_1 f d 



k=l 



N" 



N" 



1 



N" 



1 



-d< N" 



N" 



for j)i,...,p„6 J), neN. Now, using 
(3.12) 



Sj e n 



E 

p$zD : p-nmodN 



INT- 1 1 

in the form 

°p e a p (n)' 



IN- 



one deduces from (3.11) that 
(3.13) S*j£ 2 ({m e Z" | |M| < R}) 



Cf 2 meZ" 



m\\ < ||N- 1 || |I| il+(||N- 1 ||/(l- llN- 1 !!))^' 
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Thus Proposition 3.1 follows with 
(3.14) H=\neZ u 



< ( IV 



-I/MIn-I))*} 



□ 



Remark 3.2. The other method of proving Proposition ^lj is a small variation 
which gives an optimal choice of H given only D. By a theorem of Bandt Jl], |[ |ll| 
^9) cited in || (3.11)— (3.12)] there is a unique compact subset X C W such that 
X is a fixed point for the map Y 1— > {J p€D o~ p (Y), i.e., 

(3.15) X=\Ja p (X), 

p£D 

and we may take 

(3.16) H = X(1Z V =H(D). 

In some examples in Section [|, the finite subset H C Z" will be computed explicitly. 
If the representation of O^r is irreducible, an application of || Lemma 6.1] further 
shows that the finite-dimensional subspace IC(H) from ( |3.16 ) contains a unique 
minimal subspace M. with the invariance property S*M C M.. 

The following corollary is the main tool in analyzing polynomial representation. 



Corollary 3.3. Consider the 'polynomial representation of On defined by (3.7) 



and (3.S), and let H be a minimal finite subset of 1^ satisfying the properties in 



Proposition 3.1. It follows that 



(3.17) K = P (H) 

is cyclic for the representation, and thus the representation is finitely correlated. 
Defining V* € B (£) by 



(3.18) 



m£H p£D: p— n mod N 

f P (n)£fl 



a p )e CTp (n) 



/or n € -ff , i/ie commutant of the representation is isometrically order isomorphic 
to 



(3.19) 



B{K,y = IA&B{K,) 



In particular the representation is irreducible if and only if < B(/C) <T = CI. In i/izs 
case, the peripheral spectrum of a is always a finite (necessarily cyclic) subgroup 
of T, and if k is the order of this subgroup, the restriction of the representation to 
UHFjv decomposes into the direct sum of k mutually disjoint irreducible represen- 
tations. 

In general the intertwiner space between two representations of this type is given 
by p^D-(p4j). 

Proof. The identity 

(3.20) 1= V SjS*j, 



E 

/: \I\=M 



in conjunction with Proposition |3.l| , implies that all monomials e n , n € Z", are 
contained in the cyclic subspace generated by /C, and hence this space is dense in 
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L 2 (T"). Indeed, for every n G W , there is, by Proposition |3.l| , an M € N such that 
S} e € K. for all / such that |/| > M. Therefore SjS*je n e Sj (fC). An application 



of (3.20) to e„ then yields the desired cyclicity. This cyclicity is the second of the 
two properties of the subspace JC in the discussion of Section ||, i.e., ( |2.9| ). The rest 
(and some more details) follows from the discussion in Section 0. □ 

4. Classification of some polynomial representations 



If D is a given finite subset of 7L V , the set of all polynomials rrij given by (3.6) 
and satisfying the unitarity condition ( |3.3] ) and the normalization 



(4.1) m (l) = VN 

(which is necessary for the convergence of the Mallat expansion; see |^3| or [|[ 
cq. (1.37)]), forms a compact algebraic variety M.O-, and it is given as the solution 

variety of a set of quadratic equations in the coefficients an and an with n G 
D. For each point on this variety Md, the corresponding representation of Oat 



can in principle be computed from Corollary 3.3. Even the characterization of 
Md is a formidable task in general, but it has been done in the case v = 1 and 
N = 2 in g H |§ f^, ]l§ and |l| (see also g |7|)- In this section, we will 
compute the representation theory of O2 for each of the points of some of these 
varieties. We do not know if our results indicate how the generic behaviour of this 
representation theory will be, but in the examples the representations generically 
are irreducible and mutually disjoint, with exceptional behaviour on a sub-variety 
of lower dimension. 

4.1. The case with dimension v = 1. In this case, N = N 6 {2, 3, 4, . . . }. If 

rrij (z) — X^neD °n **> where m = mi for some i, unitarity ol 

(4.2) ]T \ mj (p k z)\ 2 =N, 

which is equivalent to the conditions 



(4.3) 5>«a« = l and £a«a« roJV = 

n n 

for m = 1, 2, . . . . Analogously, orthogonality of the rows in (|l.6|) leads to 



(4-4) E^^- 



for all i ^ j and a// m € Z. Finally, the normalization (4.1) leads to 
(4.5) £^ 0) =vA- 



The relations (4.3)-( |4~5| ), together with an = for n £ D, determine the 
algebraic variety M.d- Let us now restrict to W = 2, and to the case where 
the an^'s are real (this latter assumption, reality, seems conventional in wavelet 



theory). Then by (|1.14| ) 



(4.6) m 1 (z) = zf (z 2 ) m (-z), 
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where / is a monomial. By translating the father, and mother, functions by multi- 
ples of T (integral translations), we may assume that D has the form 

(4.7) D = {0,l,...,2d-1}, 

where d £ N, and with / (w) — — we have 



2d-l 



(4.8) 



(4.9) 



m (z) = ^ a k zk , 

k=0 
2d~l 



2d-l 



2d-l- 



The conditions (4.3) 
(4.10) 



a k z 
then become 



fc=0 



/J Qfe = 1, and a k a k+2r , 

k=0 k=Q 

for m = 1, . . . , d — 1 (no condition if d = 1), and 

2d-l 

(4.11) £ a fc = V2. 



0, 



(The condition (4.4) is already take n care of in (4.9).) 
In this case the maps a p in (3.10) have the form 

x — p 



.12) 



a p (x) 



for p = 0, 1, . . . , 2d— 1, and thus the solution X to the equation (3.15) is the interval 

(4.13) X = [-2d+l,G], 
and hence by ( 3.16| ) 

(4.14) ff = {-2d+l,-2d+2,...,0}, 

(4.15) K = span {e_ 2 d+i, e_ 2 d+2, ■ ■ • , e } . 



It follows from (3. IS) that the matrix for Vq relative to the basis {eo, e_i, . . . , e_2d+i} 
has the form (passing under the name "slant-Toeplitz matrix") 

(4.16) 



a 






















... 














\ 


a 2 


Oi 




a 














... 
















04 


«3 




«2 




ai 




o 




... 


















d2d- 


5 


a-2d- 


6 


02tf- 


7 


a2d- 


8 


... ai 


a 













«2d-2 


a.2d~ 


3 


a,2d- 


4 


a2rf- 


5 


d2d- 


6 


... a 3 


»2 


ai 


a 










d2d- 


1 


a,2d- 


2 


02rf- 


3 


d2d- 


4 . . . 


... a 5 


04 


a>3 


a 2 


ai 

















a2rf- 


1 


d2d- 


2 


... a 7 


06 


05 


04 


a3 



























■■■ | 0S2ei-l 


02d-2 


»2d-3 


a2d-4 


a2d- 


5 

























... 





02d-l 


a 2 d-2 


02ci- 


3 

























... 











a2rf- 


J 
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and the matrix for V* is, by (4.9), obtained by using the substitution a k 



(— 1) azd-i-k m the matrix ( 4.16 ). Note that the subspace 
(4.17) /C = span {e_ 2 d+2, £-2d+3, ■ ■ ■ , e_i} 

is also invariant under V * and Vf* , and thus under Sq and SI , but we will see in 
Section 4.1.2.2 below that this subspace is not always cyclic. 

Let us remark that the scaling relations for the father function ip corresponding 
to (4.8) and the mother function ip from (4.9) (both in I? (W.)) are as follows: 



(4.18) 
(4.19) 



1 



X?) 



a k ip (x - k) 



(-1) avd-i-W (x — k) . 



See also Remark 4.3 



Following the terminology in 
to a closer study of d < 2. 



we say that d is the genus, and we now turn 



4.1.1. The case with di mensi on v = 1, scale N — 2, and genus d — 1. I n this case , 



the second condition of ( 4.1C ) is vacuous, and the only solution of ( 4.1C ) and ( 4.11 ) 
is a = ai = 4=, so 



(4.20) 



m 0) = (1 + 2) /V2, mi (z) = (1 - z)/V2, 



which is exactly the Haar wavelet (Figure ||). The representati on sp lits in to th e 
direct sum of the two inequivalent irreducible representations in ( 1.18 ) and ( 1.19 ), 
and the restriction of each of these representations to UHF2 is still irreducible by 
|U Proposition 8.1]. This can also be checked directly: in this case, 



(4.21) 





Thus 



(4.22) 



t=0 



a b 
c d 



V* = 



a 
d 



so B (/C) CT is the *-algebra of all diagonal 2x2 matrices. Thus the representa- 
tion splits into the direct sum of two representations with the one-dimensional 
S'*-invariant subspaces Ceo and Ce_i. The corresponding maps cr on the one- 
dimensional subspaces are both equal to the identity, thus they are ergodic with 



peripheral spectrum 1, and UHF2 is dense by Corollary 3.3. Note that the states 
on 02 corresponding to eo and e_i are the Cuntz states (see (Io[ [l(], |(|) 



(4.23) 



w (S I S* J ) = 2~ 



W _i(^) = (-1)'^2-^ 



4.1.2. The case with dimension v = 1, scale N — 2, and genus d = 2. We now dis- 
play the one-parameter family (with two singular points) of mutually inequivalent 
irreducible representations of O2 mentioned in the Introduction, and we relate the 
representation-theoretic behavior to the corresponding properties of the associated 
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Values of 9 
shown: 




a: ip for 9 = \ 



1 
1 


/ 






r 


12 3 
a': V for = 7 -f 



-1 



12 

c: <p for 9 = 


3 

7 


( 




^ \ 





12 3 

c': V for 6 = ^ 




12 

b: p for 6» = i|2L (or 



b': V for = ^ (or = -f ) 




12 

d: p for 9 = 



h 
1 \ 



•\ 'J 



d': V for 9 = ^ 



2 3 

)7T 

4 



Figure 3. Father (p) and mother (?/>) functions for 9 near ^ and 
— |: Continuous cases (Case "a" = Daubechies wavelet) 
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family of wavelets on L 2 



In this case the algebraic variety defined by (4.10)- 



(4.11) is actually the circle, and may be defined by the following parametrization: 



1 



(4.24) 



a = 



a 2 = 



2V2 
1 



(1 — cos^ + sin#) , 
(1 + cos 9 — sin 9) , 



ai = — 7= (1 — cos 9 — sin ( 
2v 2 

as = — ^ (1 + cos 9 + sin f 



2V2 y " ° 2V2 

see 25 , 26 , and also 30 , 19| . Let us give a simple argument for this parametrization: 
View a — (do, a\, 0,2, 03) as a function on the cyclic group of order 4, Z 4 , and 
consider the Fourier transform on Z4: 



(») = \ L * ,mia ( m ) > a ( m ) = \ E r nmfi ( n ) ■ 

m— n— 

ual formulae 

a (m)b (m) = V] a (n)S (n) , a ( • + k) (n) = i~ nk a (n) , 

m n 

a (m + 2)a (m) = (— 1)" a (n)a (n) . 



m=0 

We have the usual formulae 



and thus 



Also a(n) — a{—n). The relations ( [4.3] ) and (4J3), together with reality of a, take 
the form 



d n d 7r 1, 



d n d n . 



and hence 



+ aia3 = a (?ti + 2)a (m) = 0, 



«(°) = 7f' 5>(»)| 2 = 1, 



a (— n) = a (n) ■<=>■ a (0) , a (2) are real and a (3) = a (1), 

£(-i) n afr)fi(n) = o. 

Introducing c = a (2) = c and 6 = a (1) we thus have c 2 + 2 |&| 2 = 1/2, c 2 — 2 \b\ = 
— 1/2, and hence c = 0, |6| = 1/2. Putting b = \e lv , the relations for a are thus 
equivalent to 



(a (0), 6 (1), o(2), a (3)) 



1 e iV 



,0,- 



,V2' 2 2 
Applying the inverse Fourier transform to this, we obtain 



1 



«2 



2V2 
1 

2V2 



1 + V2 cos 95 
1 — \/2cos<£ 



1 



«i 



03 



2N/2 
1 

2\/2 



1 + V2 sinyj 
1 — v2 sin y> 



Substituting 95 = 6> + ?f here, we obtain (4.24). 
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Returning to the representation, the operators V£ from (4.16) in this case have 
the form: 



(4.25) 



Vo = 



( ao 








°\ 


a 2 


ai 


a 








Ct3 


a 2 


Ol 


\0 








as) 



and V* 



(0.3 








° \ 




-a 2 


A3 








-ao 


ai 


-0.2 



\0 -ao/ 



If one replaces the angle variable 9 with ip, and calls the corresponding coeffi- 
cients bo, . . . , 63, and th e cor responding matrices Wq, W£, the corresponding map 
p: M 4 -> M 4 given by ( |2. 14| ) , 



(4.26) 



1 

£ 

z=0 



is defined by a 16 x 16 matrix relative to the basis 
(4.27) e ,o, e -1, . . . , e 0: _ 3 , e_i, Q , e_i _i, . 

of M4. This 16 x 16 matrix has the form 

/A A 2 \ 
r42Sl Ui A 3 \ 

(4 - 28) Uq A 2 \ ' 

\0 Ax A 3 J 

where the 4x4 matrices Ai are given by 

Ao = boV + hVi, A x = bxVo - b 2 V x , 
A 2 = b 2 V Q + bxVx, A 3 = b 3 V Q - b V x . 



e-3,-3 



(4.29) 



Thus one can compute the eigenvalues of p by computing the eigenvalues of the 
matrices A , A 3 , and ( ^ ^ 3 ) • H <p = 6 the result is (we call p = cr in this case as 
usual) 

s 9 cos 6 1 



(4.30) 



Eigenvalue of <r 



Multiplicity 



1 



sin ( 



— sm( 



1 



A} is 



Hence, the dimension of the eigenspace {A | cr (A) 

( 3 if 9 = f 

(4.31) <2 if0=3f 

I 1 otherwise. 

These numbers are then the dimensions of the commutants of the corresponding 
representations. Since the only C* -algebras of dimensions 1, 2, 3 are C, C 2 , C 3 , 
it follows that the representation of 2 splits into 2 inequivalent irreducible repre- 
sentations if 9 = 4^, into 3 inequivalent irreducible representations when 9 = 
and the representation is irreducible for all other 9. We note that the peripheral 
spectrum of cr is nontrivial only if 9 = ^, when —1 is an eigenvalue of multiplicity 
1. Thus the representations for generic 9 ^ {f : a l so have irreducible restric- 
tion to UHF2. Finally, if one considers the case 9 ^ <p, one can compute that 1 
is an eigenvalue for p if and only if {9, ip} = { § , 4 f }, a nd the dimension of the 
corresponding eigenspace is then 2. We recall from ( |2.14 ) that solutions A 7^ to 
p (A) = A correspond by lifting to operators on L 2 (T) which intertwine the two 
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associated 02-representations and for 9 and tp, respectively. Hence the 
representations for generic points 9 £ {§,%-} on the circle are all mutually disjoint 
by ( 2.13 )-( 2~15|), b ut if {9,ip} = {§, 4^}, the intertwiner space is 2-dimensional. 
See Section 4.1.2.3 for more details on the latter. 

A second immediate observation on ( |4.24| ) is that at the four points 9 — 0, -|, n, 
we have two of the four coefficients vanishing with different pairs in the four differ- 
ent cases, so those four cases are closely connected to four modified Haar wavelets, 
illustrated in Figures |andf§ A more subtle fact, to be described below, is that 
it is only the two cases 9 = 5 and 9 — 4^ on the symmetry axis where the cor- 
responding 02-representation on L 2 (T) fails to be irreducible. The case 9 — -| is 
degenerate in a sense illustrated in Figure [|. We will relate the resulting degenerate 
decomposition at 9 — ^ of the subalgebra UHF 2 C O2 to the wavelet properties. 

Let us now consider the two exceptional points 9 = \ and 6 = separately. 

4.1.2.1. The case 9 = 4f. When 9 = a = a 3 = 0, ai = a 2 = 1/V2, so 

(4.32) m (z) = (z + z 2 ) /y/2, ip {x/2) = <p (x - 1) + <p (x - 2) , 

(4.33) mi (z)= (-z + z 2 ) /V2, i/;(x/2)=-tp(x-l)+ip(x-2), 

with the scaling relations indicated for the father function tp, and the mother func- 
tion tfj, respectively; see Figure |[ 

This is a simple transform of the Haar wav elet (Figure |2|) , and the representation 
theory becomes similar: defining Si by (F9) and transforming the representation 
by 73 (ill) eU(2),i.e, 

(4.34) T = (S - Si) A/2, T x = (S + Si) /y/2, 
we obtain 

(4.35) T i{z) = 74{z 2 ) 1 T 1 (,(z) = z 2 i{z 2 ). 

By the computation in eqs. (8.1)-(8.2)], if U is the unitary operator given by 
multiplication by z , then 

(4.36) U*T UZ(z)=t;(z 2 ), U*T 1 U(i(z)=zti(z 2 ). 

By EL Proposition 8.1], L 2 (T) splits into the two irreducible subspaces spanned 



by{l, 



z, z 



} and {x 



}. Applying U to these, we obtain the two irre- 



,z 



2 , z 3 } commutes with the 



The respective products of P by these 



ducible invariant subspaces corresponding to the original representation 

(4.37) span {z , 1, z, z 2 , . . . } and span {z~ 2 , z~ 3 , . . . } 

(overbar for closure). We see that the projection P onto the overlapping four- 
dimensional S*-invariant subspace JC — span |1, z 
projection onto the first two subspaces. 
projections are 



(4.38) 



and these two matrices span exactly the eigenspace of cr corresponding to eigenvalue 
1. Also, each of the two subrepresentations has irreducible restriction to UHF 2 , 
confirming the fact that the peripheral spectrum of u consists of 1 alone. 



(1 












(0 








°\ 





1 








and 
































1 





V° 








<v 












V 
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4.1.2.2. The case = §. When 9 = f , 

(4.39) a Q = a 3 = 1/V5, d = a 2 = 0, 
so the associated low/high-pass filters and scaling relations are: 

(4.40) m (z) = (1 + z 3 ) /V2, <p (x/2) = <p (x) + <p {x - 3) , 

(4.41) mi (z) = (1 - z 3 ) /V2, ^) (x/2) = (a;) - (a; - 3) . 
See Figure [| for the graphs of the corresponding y> and '0. Applying the unitary 



€ U (2) to this representation, we transform it into the representation 

3 



i/ii 
TSu-i, 

with mo(z) = 1, mi (z) = z 3 . We have already noted in (4.31) that the fixed 
point set of er is three-dimensional in this case, and indeed, by |t|, Proposition 8.2], 
this representation decomposes into 3 mutually disjoint irreducible representations 
given by restriction to the 3 subspaces 



(4.42) 



span { j 



3n 



= 0,1,2,...}, span{ 
span {z fc | k not divisible by 3} 



1. 



•} 



The restriction to UHF2 is still irreducible on the first two subspaces, while it 
decomposes into the two irreducible subrepresentations on 



(4.43) 



span { 



y 3fc+l 



k £ Z} , span { 



3k+2 



kez} 



on the third subspace. Again the projection onto each of these subspaces commutes 
with P, and hence the eigenspace of er corresponding to eigenvalue 1 is spanned by 
the three projections 



(4.44) 



respectively, confirming that the eigenval ue 1 has multiplicity 3 in this case. Fur- 
thermore, if U is the unitary operator ( 2.21 ) on span {z k \ k not divisible by 3} 
that implements the gauge automorphism r_i there, we have UTi = —TiU. Hence 

(4.45) U (z 2 )) - - (UO (z 2 ) and U (z 3 £ (z 2 )) = -z 3 ([/£) (z 2 ) 

if £ is in this subspace. This unitary U from ( |2.21 ) has to fix the two subspaces 
span { z 3k+1 I k 6 Z} and span {z 3fe+2 | k G Z}, and U 2 = 1, hence it is clear that 



























(0 








°\ 




























, and 





1 








































1 





\o 












{0 




















0/ 



(4.46) 



PUP = ± 



/o 











It is easily verified directly that this is the eigenvector of er corresponding to eigen- 
value — 1. This means that the group from (2.18)-( f2.20| ) in this case is Z 2 , if JC is 
taken to be span{e_i, e-2} and ip the trace state on B (/C). 

In conclusion we note that this 02-representation tv^ 8 \ 9 — ^, as well as its 
restriction to UHF2 has a decomposition into irreducibles which sets it apart from 
the other representations when 9 ^ -|. We will see in the beginning of Section 
4.1. 2. 5| that if 9 — 0, w, or then the wavelet is still of Haar type, i.e., ip is of the 
form tp = xi where / is an interval of unit length. The position of the interval I 
varies (see Figure [|) in the three cases 9 = 0, 7r, or while the mother function 
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is common for two of them, 9 = 0, 7r, and = — = — ^W. All 

three satisfy ^> (3 — x) — —ip (x). But, if 9 — ^, then the nature of <p is somewhat 
different. From ( |4.40D , we see that ip (§) = ip (x) + cp (x — 3) ; and by [jl||, then ip 
must have the form ip = where S is a compact subset C [0, 3] with non-empty 
interior. It is determined by the identity 

(4.47) 2S = SU(S + 3) 

(see Figure ||). It follows that S = [0,3]. In fact, iteration of (4.47) leads to 
the following representation which characterizes points x in S: x = ^2k = idk/'2 k , 
dk = 3£fe, € {0,1}. Hence, using base 2 for the unit interval [0,1], we get 
S = [0,3]. The derivation of from the first Haar wavelet p^ — X[o,i] ls a 

special case of the substitution 



(4.48) 



to (z) 



m 



(z 3 ) , or generally, mo (z) 



toq 



(z 2 p +1 ) 



If too is an arbitrary low-pass filter with scaling function ip, then the argument 

--p((2p- 
Hence 



from Remark 4.2 shows that pi p +\ (to) 



will determine the scaling 



function for the substitution m (z 2p+1 ^ 
(4.49) 

(f) 



2p+ l 1 ^ V 2 P + 
In our circular family, we have Toq 



\/2p+T 



Ml 



Section 4.1.2.4 and Remark 4.3 



» 



(z 3 ). See further discussion in 



4.1.2.3. Intertwining of the cases 9 = ^ an d 9 = 



tion of these two representations. By (3.7)- 



(4.50) S > e TI 



1 ( \ 

^= ( e l+2n + e2+2n) ■ 



Let us summarize the descrip- 
8|) we have 

1 , . 

^= (-ei+2n + e2+2n) , 



span{e_2,e_ 3 , . ..} , 
(e 2 )/v / 2, 



(4.53) 



and the irreducible invariant subspaces are 

(4.51) H + 2 = span{e_i, e , ei, . . .} , Hj 1 = 
Similarly 

(4.52) S§e n = (e 2n + e 3+2n ) /V2, sfe n = 
and the associated three irreducible invariant subspaces are 

HI = span{e , e 3 , e 6 , . . . } , HI = span{e_ 3 , e_ 6 , e_ 9 , . . . } , 

J = span{. . . , e_ 4 , e_ 2 , e_i, ei,e 2 , e 4) . . . } . 

We have noted that UHF 2 is not weakly dense in the last representation but 
it is so in the first four, so the last representation cannot be equivalent to any of 
the former four. Also the representation on H e + is disjoint from that on H_ by 
(H Theorem 2.7], for 9 = ^ and for = ^. So the remaining possibility is that 

3tt tt 

the representation on 7i ± 2 is unitarily equivalent to that on H± . Inspection of the 

expressions for Sfe n makes it plausible that the representation on H+ is equivalent 

to that on HI , and that that on H2 is equivalent to that on H\ , and indeed, if 
one defines an isometry U by 

(4.54) Ue n = e- 3n - 6 
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-1 



1 2 

ip for e = 



-1 

3 



1 2 

iyS for 9 = 7r 



-l 
3 



V? for = 12- 
(cf. (p|)) 



1 1.52 

ip for 6» = 

(cf. (HD) 



-1 

3 



1 1.52 

ij) for 9 = 7r 
(cf. flOp )) 



-1 

3 



1 1.5 

i/j for = 

(cf. @) 



3tt 
2 



FIGURE 4. Father and mother (-0) functions for 
multiples of The symmetry ?/>(3 — a;) = —ip(x) 



equal to 



then £7 1 3» from W^ 2 to and (7| 4^ from H.1 to arc unitary operators, 
and one computes 
(4.55) 
(4.56) 



US^ e r 



°-6n-9 + e -6n-12j 
-e -6n-9 + e -6n-12) 



Hence £7 intertwines the two representations, and if U is restricted to 7i ± 2 one 
obtains the expected unitary intertwiners 

(4.57) U X :HJ — >H f , U 2 : Hj — >H%. 

Now Uie-i — e_3, U2&-2 — eo, and hence 

/0 y 0\ 
(4-58) P(xU 1+ yU2)P= q q q q 

\0 1 0/ 

for a:, y e C, and this is exactly the fixed point set for the map 

1 

(4.59) M 4 9 4^^W' 

i=0 



as it should be by (2.13)-(2.1E). By |6j, Theorem 5.1], these solutions A correspond 
to operators which intertwine the two representations. 
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12 3 

First approximation to ip 



12 3 

Second approximation 



12 3 

Third approximation 



12 3 

Fourth approximation 



1 

Weak L 2 limit: 



3 



(cf. (MoD ) 



^ for e = § 



1 2 

ip for = § 

(cf. dub) 



Figure 5. Father (ip) and mother (ip) functions for 9 = with 
cascade-algorithm approx imations of father function y>: See dis- 
cussion in Section |4.1.2.6| 



4.1.2.4. Additional remarks on singular points and cycles. Recall from Jl^, Theorem 
6.3.6] and Theorem 3.3.6] that in order that ipj.k (x) — 2~iip (2~^x — k) shall 
be an orthonormal basis for L 2 (R) and not merely a tight frame, it is necessary 
and sufficient that the set 

(4.60) |z e T | \m (z)\ = v 7 ^} = {z e T \ m a (-z) = 0} 

does not contain a nontrivial cycle for the doubling map z i— > z 2 , i.e., a finite 
cyclic subset unequal to {1} invariant under the map z i— » z 2 . Inspection of the 
polynomial (z) in (4.8) in the present case ( 4.24| ) reveals that the condition 
above is fulfilled for all 6 S T with the sole exception 



(4.61) 



6 = tt/2, 



where TOq (z) is given by (4.40), and thus the set ( 4. 60] ) consists of the three cube 



roots of 1. Indeed, the presence of a nontrivial cycle on T under z 



imply, by ( |4.60| ) and the fact that is a third-degree polynomial, that m^' ( ■ ) 
is a scalar multiple of (z + 1) (^z + e l ^~^j (^z + e^T 1 ^ = z 3 + 1, and this is precisely 



would 

,(*) 



the case 9 = \ in (4.61). (See details of the argument below.) It is interesting that 



this is the case of Section 4.1.2.2 where the decomposition theory of the associated 



representation is most singular, being the only case where the restriction of one of 
the subrepresentations to UHF 2 has a nontrivial cyclic structure. 
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Let us give a more detailed justification of the statement above. First note that 
cycles on T are not subgroups of T but rather cyclic orbits on T under the z i— » z 2 
action of one of the cyclic groups k = 1,2, ... . Such a cyclic orbit with k 
distinct points z\,...,Zk must be of the form z\ —* Z2 —* ■ • • — > Zk — * z\ , where 
z i+l = z 2 if % = 1, 2, . . . , k — 1, and z\ — z x . Hence points c in an orbit C'k must 

satisfy c 2 = c, and each c must be a (2 fc — l)'th root of 1. Different orbits must be 
disjoint, and their union will be invariant under z^tz 2 acting on T. The converse 
is not true. For example, the subset {1,-1} C T is invariant under z i— » z 2 while 
not a cycle, and not even the union of cycles. Note also that we can have different 
(2 h — lj'th roots c of 1 defining different cyclic orbits for the same fc. If fc = 1 
or fc = 2, then in each case there is only one orbit, but if k = 3, there are two 
choices. Since m{, 9 ' for each 9 is a polynomial of degree 3, the cardinality of a cycle 
contained in the set ( 4.60| ) is at most 3. Thus, if z is contained in such a cycle, 



we must have one of the possibilities z = z, z = z, z = z. Hence the cycles of 
length at most 3 are {1}, {kj,a; 2 } where ui = e z ~^ , {p,p 2 ,p 4 } where p = e l ^~ , and 
{p,p 2 ,p A } — {p 6 ,p 5 ,p 3 }. But as mo (— 1) = always, (z + 1) is always a factor 
of mo (z), and since the cycle should be different from the trivial cycle {1}, we 
are reduced to the case {a-i,oj 2 }. The other cycles would make m^ divisible by 
a polynomial of degree at least 4, which of course is impossible. Thus we are left 
with the case 

(4.62) m (z) = ± J] (c fc + z) = ± (l + z 3 ) , 

which is exactly the case 9 = | . 

Note, more generally, that the wavelets which arise from substitutions, as defined 
in (4.48)-(4.49) with filter function m^ (z) = m (z 2p+1 ), will have those addi- 



tional cycles Cfc which are contained in the (2p + l)'th roots of 1, {z e T | z 2p+1 = l}. 
We will show in a forthcoming paper that this leads to a decomposition of the rep- 
resentation of O2 associated to mjj over the new cycles. 

It is interesting to note that the same cycles as described above arise in a differ- 
ent context in Q in connection with a family of discrete series of representations 
of On- These representations are called permutative representations, and the cy- 
cles represent the finite decompositions of irreducible representations of On when 
restricted to UHFat. 

4.1.2.5. The symmetry 9 tt — 9. Note that the two points 9 = and 9 = tt are 
interesting in that the representation theory is regular, but these points correspond 
to mother and father functions which are simple rescalings of those of the Haar 
wavelet (see Figure |i|) : 

(4.63) 

'mo(z) =(z 2 + z 3 )/V2 6) = 7r fm (z) = (1 + z) /y/2 
mi(z) =(l-z)/V2 \mi(z) = (z 2 - z 3 ) /y/2 

Thus the representation of O2 is very sensitive to simple rescaling of (p and ip. In 
fact the mother function ip is the same in the two cases 9 = and 9 = tt, and this 
common ip has the following symmetry property tp (3 — x) = (x), which in turn 



= 0: 



is a special case of a more general reflection symmetry (4.67) to be discussed in 
Proposition 4.1(a|) below. 
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The symmetry h- > ~n is a special case of a symmetry 9 i— > -k — 9, which we will 
now analyze further. If this transformation is substituted in (4.24), we note that it 
corresponds to the following reversal: 

(4.64) (a ,a 2 , a 2 ,a 3 ) i — > (a 3 , a 2 , ai,a ) ; 
or equivalently, 

(4.65) m^~ e) (z) = rn£ ] (z^ 1 ) z 3 . 

The following proposition shows that the 9 i— > 7r — reflection applied to to 
implements the x h- > 3 — x transformation on the scaling function ip (see Figure 
[|a,b). It is interesting to note that, despite this left-right mirror symmetry of the 
graphs in the family of scaling functions p^ s \ the two associated representations 
of 02 on L 2 (T) which correspond, respectively, to 9 and ir — 9, are not unitarily 
equivalent, by the results above, except of course at the two fixed points ? and ^ 
for 9 i — > 7r — where the representation theory also happens to be exceptional. See 
subsections 4.1.2.1 and 4.1.2.2 above. 

Proposition 4.1. Let m (z) be the filter functions indexed by 9 and correspond- 
ing to the given coefficients in the family (4.24). Let cp( e > (x) be the associated scal- 
ing function (alias, father function) and the mother function corresponding to 
the pair (m^\ of low/high-pass wavelet filters. Let sf^ be the corresponding 

operators from (3.4). 

(a) The symmetry relations 

(4.66) ip {7r - 0) (x) = ip {e) (3 - x) , 9 € [-7T, tt] , x e R, 

(4.67) V ( "~ e) (a:) = -^ {6) (3 - x) , 
are mZid. 

(b) The corresponding representations tz^ and tt^^ 9 ^ (given by tv^ (sj = ) 
satisfy 

(4.68) IT7T W = ( 7r (7r ~ 9) O Tfl o\)W, 



.0-1. 



where (Wf) (z) = z f (z ), awe? 77 1 o A * s ^ e automorphism of 2 given 
in^)forg=(l\). [o 

Proof. Introducing z = e~ lu) , w 6 R, the identity ( 4.65] ) above reads 

(4.69) m^~ S) (lu) = e l3w m^ (-w) = e^m^ (lu), wel. 

Generally for third degree, the correspondence too «-> <p is given by the following 
functional identity in L 2 (R): 

(4.70) if (x/2) / V2 — a Q ip (x) + a\ip (x — 1) + a 2 ip (x — 2) + a 3 <^ (x — 3) , 

and the boundary conditions, tp(0) = Lp (3) = 0, i.e., </? is uniquely determined by 
these conditions and the normalization (p (0) = (2ir) 1 . See |27j for details. This 
applies to both the pair [m^\ (f^^j and the pair [rri^^ e \ ^( 7r_0 )J , so we get 

(4.71) p^ (x/2) A/2 

= a 3 ^ (,r_e) (a:) + a 2 <^ ( ^ 9) (a - 1) + ai<^- 9) (x - 2) + (a - 3) . 
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As noted, ( • ) is the unique normalized L 2 ( 



ject to ip^~ 6 ) (0) 
oW (3 



(3) 



0. 



to rn,Q J , then a direct substitution x 



x i— > ip y 
(4.72) 



a;) satisfies (4.71 



solution to this identity, sub- 
But, if y>W is the solution corresponding 
-> 6 — a; shows that the mirrored function 
and we conclude from the uniqueness that 



(x) = ^ (3 - 



x) , x € 



as claimed in the Proposition. The proof of (4.67| ) is similar, or see Remark [f.2 
below. We resume the proof of Proposition 4.1 



(bh after the following remark. 



Remark 4.2. Proposition |4.l| (|a|) may alternatively be proved from the Mallat al- 
gorithm as follows: If f^, -0^' are the father and m other functi ons a t the angle 
6, and the transformation 9 i— > 7t — 6 is used on ( [4.24 ) , we obtain ( 4.64 ) and ( 4.65 ) 
as before, i.e., 



(4.73) too (z) i — ► z 3 too (z) = mi (—z) and mi (z) i — ► too (—2) . 

Applying the Mallat algorithm <p{t) = (2tt)~5 (m (e~ rf2 ~ fc ) /v^J, we ob- 

tain 

(4.74) <^ r? ) (i) = e^ t3t (^W (-*) , 
and thus by Fourier transform, 

(4.75) V {7! - B) (x) = y>W (3 - x) , 



which is (4.66). On the other hand, 



(4.76) 



and so 



(x) = (-!) fc a 3-fc^ (2a: - fc) , 



(4.77) (a:) = (-1)" a^V^^ (2s - k) 

k 

= V^(-l) fe 4V 9) (3-(2z-fc)) 
fe 

= v/2 £ (-l) 3 ~ fc a< V< 9 > (2 (3 - x) - k) = (3-x). 



which is (|4.67| ) . 

It is important to note that the infinite product argument works even if tp^ [x] 
is not continuous in x. Since 



(4.78) 



2, w e 



to (uo) + m (lj + 7r) 

it is known that the infinite products 

00 00 

(4.79) (a*)"* n (£) , (Jbr)"* 2"M 9) (f ) II ( 



fe=i 



fc=2 



are well defined and represent (p( e \ tp( e \ where <p( >,ip( > e L 2 (R) 12 
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Proof of Proposition [i.l| (|b|) . Let us consider the two operators Sq^ and Sq ^ in 
L 2 (T) individually, and as part of a pair of (^-representations. While the two O2- 
representations are inequivalent, the two So-operators alone are unitarily equivalent. 



This follows from the general fact that any operator of the form (L9) coming 
from a wavelet is unitarily equivalent to the shift of infinite multiplicity by [g| 
Lemma 9.3]. The explicit intertwiner can also be calculated directly as follows: Let 
m (z) = ao + a\Z + • • • + ar>z D , m' (z) := z D m (z^ 1 ) , and define three operators S, 
S', and W (acting on f e L 2 (T)) by 

Sf(z) :=m(z)/(z 2 ), S' f (z) := m' (z) f (z 2 ) , Wf (z) := Z -° f (z" 1 ) . 

Then W : L 2 (T) — » L 2 (T) is a unitary intertwining operator for S and S", i.e., 

WS = S'W 

holds, as can be verified by a direct calculation. The minus sign in the second 



symmetry formula (4.67) is still reflected in the (^-representations as follows. Let 
D = 3 and m — mj 9 ' , and consider the two 02-representations , n^~ ^ , i = 0, 1. 
We then have m'f^ 8 ^ (z) — —z 3 m^ and thus 

(4.80) WS { Q e) = S^~ e) W, WS{ 9) = s[*~ $) W, 

where again = (sj. Hence W intertwines the ^-representation with 
the (it — ^-representation Tv^^ e \ modified by the automorphism of 2 induced by 



9={h°i) eU(2);see(£3). □ 



4.1.2.6. Contin uity o f scaling functions: Stability interval. Historically the special 
case 8 — in ( 4.24 ) was discovered first. In that case, 



l + V^ 3 + ^/3 3-v^ 1-V3 

a °~^7T' ai ~^7T' a2 -^7T' a3 "^vT' 

which is the (by now) well known Daubcchics wavelet; sec Figure |a,b JH, Chapter 
6] . It was analyzed further in |27| , where it was shown to have scaling function ip ( • ) 
continuous and one-sided differentiable in x, support on [0,3], ip (0) = (p (3) = 0. 
It is left-differentiable at every dyadic x, but it is not right-differentiable at any 
dyadic x in [0,3). Since Daubechies established continuity by a matrix spectral 
estimate (see jl2[ Theorem 7.2.1]), it follows from her estimates that the scaling 
function ip^ (x) will also be continuous in an open interval containing 9 = ^ . It is 
interesting to note that Daubcchies's spectral estimation involves the two matrices 



V* , i — 0, 1, given in ( 4.25 ) above. The discussion in our previous section indicates 
that the stability interval in the 9 variable, 9q < 9 < 9\, must have 7r < 80 and 

2 



h < 



The pictures of the scaling function in this paper are generated with the aid of 
the cascade algorithm described in JlJ], Section 6.5]. 

For uniform convergence of the cascade approximants to ip, one has to assume 
that cp is Holder continuous jl2|, Proposition 6.5.2]. Figure || shows clearly that 
this uniform convergence may fail abysmally even when ip is a simple step function. 
However, we see from Figure |^ that the cascade approximants converge in the 
distribution sense, and even in the weak-L 2 sense, to (p when 9 = \- 

The other assumption in Daubechies 's cascade approximation is the orthogonal- 
ity of Z-translates, in the form |l^, (6.5.4)-(6.5.5), p. 204], and, as we will discuss 
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in Remark 
by Section 



| 



below, that fails when 9 



but is satisfied at all other values of 9 



2.4 above. 

More importantly, Daubechies states in |l2], Chapter 6 footnote 9 and Section 
6.3] that, even if <p is not assumed continuous, we still have L 2 (R) norm convergence 
of the cascade- algorithm approximation, as long as the Z-translates are mutually 
orthogonal; and, as we noted, this on 
be discussed in a forthcoming paper 



orthogonal; and, as we noted, this orthogonality holds whenever 9 ^ ^. This will 



Remark 4.3. Since ip^ 2 ) = -|x[o,3]i it is geometrically clear that the Z-translates 
of ip^) in L 2 (R) will not be mutually orthogonal (see Figure ||), and we have 
shown in Section 4.1.2.4 that ipi^) is the unique scaling function in the family 



|<p( e )} which does not have orthogonal Z-translates. The cascade algorithm, which 
is used in generating the present graphics, is based on an iteration of ( 4.18] ) but is 
also closely connected to iteration of F * in ( 1.1 2j ) . Let 



(4.81) 



Ck 



%(f) 



dx. 



In t he ca se when {(p ( • — k)} k€Z is an orthonormal basis, we get Ck — ak, k G Z, 
by ( 4.18 ); but, in general, we have a discrepancy Ck ^ ak which leads to a rather 
poor approximation with a^-cascades. For a more explicit estimate we need the 
following: 

Lemma 4.4. Let mo be a low-pass wavelet filter with corresponding scaling func- 
tion (p and suppose that the It-translates of ip are orthogonal. Let Lp p be the scaling 
function corresponding to the substitution mo (z 2p+1 ), and let 



V2 



<r,> 



(x - k)ip p (j 



dx. 



2 1 

< 



2p+ 1 



Then 

(4.82) £|c?> 

k 

Proof. From Uj or §, Proposition 12.4], we have 

(4.83) ]T|<^ + 27rO| 2 ^. 



/el 



Since tp p (uj) — (p ((2p + 1) uj), we conclude that 

(4.84) ^((2p+l)(c + 27r0)| 2 <^. 
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Figure 6. Correlation coefficients for 9 = ^: The correlation 
coefficients Ck, for which J2k=-2 c 1 = If < §' as com P are d to 
the scaling coefficients a^, for which X\=o a fe = i - 



This second summation is just one of the 2p + 1 residue classes for the full Z 



summation in (4.83). But the formula for Ck yields 



/ e ifc "m ((2p + 1) w) V |0 ((2p + 1) (w + 2tt/))| 2 
= 2tt^ \m ({2p+l)u)\ 2 (^\tp{{2p+l)(Lu + 2Trl))\ 2 ^ dw 

i r 27r 2p / \ 2 

(2p+l)27T J / l TO oH| 2 E^27r|0(a ) + (j + (2p+l)O27r)| 2 J du, 
~ f2 D +n27r / l m ° ^)| 2 EE 2ff l^ + C? + (2P + 1) 2tt)| 2 rfcj 

I P T ) JO =0 ; z v v ' 

n 

1 /" 27r 



2p 



1 



f 



2p + 1 2p + 1 



□ 



We now illustrate this for (^( 2 ). Sin ce y( 2 ) = ^X[o,3]> it is eas Y to compute 
exactly the correlation coefficients Ck of (|4.8l[) . The nonzero coefficients are: 

c-2 = c 5 = 1/9V2, c-i =c 4 = 2/9V2, and c = a = c 2 = c 3 = 1/3V2, 

which should be compared with ( 1.3Sp . They are also illustrated in Figure ^, and 



a comparison with Figure g suggests that replacing the a^s in the cascades with 
the Cfc's might possibly lead to a better approximation. Good approximations are 
not known in the non-orthogonal case. For more details, see |l2|, pp. 204-206]. 
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The Cfc numbers are those which may be inserted into the spline approximation 
that is also discussed in |l2|, pp. 206-207] to build in tight frame parameters in the 
approximation. 

The problem with this substitution of the Cfc's into the cascade algorithm is that, 
in the non-orthogonal case, we will have (see Lemma [O]) J^k \°k\ 2 < 1- Compare 
this to the normalization property J^k \ a k\ 2 — 1 from ( |4.3| ), or (4.10) in the special 
case. 



5. Conclusions 

We have demonstrated how a representation-theoretic approach to the construc- 
tion of compactly supported wavelets in M. d leads to: 

(i) a coordinate-free display of the examples, 

(ii) a finite-dimensional matrix algorithm for computing irrcducibility properties, 

(iii) a formula for decomposition into orthogonal sums of irreducibles. 

The theory is illustrated in the simplest cases where the power of the representation- 
theoretic approach comes into play. 
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